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^ ■ We discuss the exclusive radiative decays B — > K*^, B pj, and B — > uj'y in QCD 

factorization within the Standard Model. The analysis is based on the heavy-quark 
limit of QCD. Our results for these decays are complete to next-to-leading order 
in QCD and to leading order in the heavy-quark limit. Special emphasis is placed 
on constraining the CKM- unitarity triangle from these observables. We propose 
a theoretically clean method to determine CKM parameters from the ratio of the 
B — > plv decay spectrum to the branching fraction of B ^ pj. The method is 
based on the cancellation of soft hadronic form factors in the large energy limit, 
which occurs in a suitable region of phase space. The ratio of the i3 — > p7 and B — > 
i^*7 branching fractions determines the side Rt of the standard unitarity triangle 
with reduced hadronic uncertainties. The recent Babar bound on B{B^ — > ^"^7) 
implies Rt < 0.81 (^/1.3), with the limiting uncertainty coming only from the SU(3) 
breaking form factor ratio ^. This constraint is already getting competitive with 
the constraint from Bg-Bg mixing. Phenomenological implications from isospin- 
breaking effects are briefly discussed. 



1 Introduction 

The radiative transition 6 — > S7 is one of the most important processes for the study of 
flavour physics. As a flavour- changing neutral current interaction it is a genuine quantum 
effect within the Standard Model (SM) and has a high sensitivity to new dynamics at 
short- distance scales. The cleanest way to probe 6 — S7 is the measurement of the 
inclusive decay B Xg'j, where the impact of strong interactions is well under control 
(see PP for a recent review). For exclusive channels such as B K*'~f, which depend on 
hadronic quantities describing the hadronization of the final state quarks into a single 
K*, a theoretical treatment is more difficult. At present, the decay B —>■ Xgj already 
yields strong tests of the SM and valuable constraints on its possible extensions. 

In contrast, not much is currently known experimentally about b —>■ d'-f transitions, 
the Cabibbo-suppressed counterparts of b ^ 57. They depend on the less well deter- 
mined weak mixing parameter Vtd, rather than Vts, and could be differently affected 
by new physics. For these reasons a measurement of 6 — ^7 will be very important. 
However, the inclusive measurement of B X^'y, theoretically prefered, appears almost 
impossible because of the dominating background from B Xs7. Therefore, exclusive 
channels such as B ^ fr/ and B wy become the only way to access b ^ d'y transitions 
in the foreseeable future. 

The CP-averaged branching ratios of exclusive radiative channels are measured to be 



5(5° ^ K*°-f) = (4.01 ± 0.20) ■ 10"^ (1) 
5(5+ ^ ir*+7) = (4.03 ± 0.26) ■ 10"^ (2) 

and bounded with 90% confidence level by Babar as jS] 



5(5° - 


cj°7) < 1.0 • 


■ lO"'' 


(3) 


5(5° - 


p°7) < 0.4 • 


• 10-^ 


(4) 


5(5+ - 


p+7) < 1.8 • 


■ 10"^ 


(5) 



The corresponding results from Belle read ^ 



5(5° - 


-.^°7) 


< 0.8 ■ 


■ 10"^ 


(6) 


5(5° - 


-A) 


< 0.8 ' 


■ 10"^ 


(7) 


5(5+ - 


.p+7) 


< 2.2 


■ 10-^ 


(8) 



Even though a theoretical treatment of the exclusive decays 5 p'j, 5 uj'y, and 
5 — s> K*'y is more challenging than of the inclusive modes, there are circumstances that 
help us to make this task tractable and that will eventually yield useful phenomenolog- 
ical results. First, recent studies of exclusive hadronic modes in the heavy-quark limit 
have led to a better understanding of the strong dynamics of these decays |5|5|7j by 
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establishing factorization formulas in QCD. For the decays B V'y |6|7j this approach 
resulted in particular in a calculation of light-quark loop amplitudes that before consti- 
tuted an uncontrollable source of uncertainty. In addition it became possible to extend 
the computation of -B ^ V'j amplitudes systematically to next-to-leading order (NLO) 
in QCD IRITIS! , improving on previous analyses Second, the impact of hadronic 

form factors, which dominates theoretical uncertainties, can be reduced by taking the 
ratio B{B p^)/B{B K*l)- The ratio of the corresponding form factors is equal to 
unity in the limit of SU (3)-flavour symmetry and the hadronic uncertainty is reduced to 
the effect of SU{'i) breaking, which still needs to be estimated. Furthermore, the ratio 
of the B fTf and B K*'~f branching fractions is, at leading order in as, directly 
proportional to the side Rt in the standard unitarity triangle (UT), where 



Vt 



ts 



(9) 



Here A, p, and f] are Wolfenstein parameters. Having the complete NLO result for the 
decay amplitudes in i? — > V'-y at hand, we can calculate corrections to their relation 
with Rt and evaluate the implications in the {p,ff) plane |ll|12|13j . 

Another possibility to reduce hadronic uncertainties consists in taking the ratio of 
B plv and B ^ p'y branching fractions. Using relations between the B ^ p form 
factors in the large energy linit, it can be shown that this ratio is free of long-distance 
QCD effects in a certain region of B plv phase space. The form factors cancel in 
this situation, up to calculable 0{as) corrections, which leads to a model- independent 
relationship of B pi u and B ^ p'y observables to the CKM quantity 



VudVub 



VuVt, 



2 -2 I -2 

- (10) 



(1 - -pf + r/2 

It is the purpose of this paper to investigate how B V'y decays can be used 
to constrain the parameters of the unitarity triangle. Such constraints simultaneously 
provide a test for new physics. The various sources of uncertainty will be discussed 
in detail in order to quantify the potential of these important decays. In section 2 we 
recall the analysis of B —* V'j decays at next-to-leading order within the framework of 
factorization in the heavy-quark limit. The extraction of CKM parameters based on the 
ratios B{B — > pj)/B{B K*'y) is the subject of section 3. In section 4 we discuss how 
theoretically clean information on CKM quantities can be obtained from combining a 
measurement of B{B — ^ p^'j) with a Dalitz-plot analysis of B —>■ plv decays. Section 
5 contains an update on observables of isospin breaking in i? ^ V'^ and section 6 is 
devoted to a discussion of the decay mode B uj^. We present our conclusions in 
section 7. 



2 B ^ ^7 at NLO in QCD 

Let us briefly summarize the basic formulas relevant for the analysis ofB^ V'-f at next- 
to-leading order in QCD. For more details we refer the reader to |(j|14j . The effective 
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weak Hamiltonian for b ^ s'y transitions is 

2 8 

^^// = % E ( E c^Q'. + E c^Q^) (11) 

p=u,c ^ j=l i=3 ' 

where 

A^.^) = V;,V,, (12) 

The relevant operators are the current-current operators Q12, the QCD-penguin oper- 
ators Q3,,,6, and the electro- and chromomagnetic penguin operators Qt^s- The most 
important contributions come from and Qj^sy which read 

Ql = {sp)v^A{ph)v-A Ql = {SiPj)v-A{pjbi)v~A (13) 

Q7 = s.a^-'il + 75)6. F,, Qs = /^m, s,a'^^(l + ^,)^bj G% (14) 

The impact of penguin operators Q3, ■ ■ - Qq is very small for most applications, but will 
be included in the numerical results presented below. The effective Hamiltonian for 
b ^ d'y is obtained from ()11|) - ()14|) by the replacement s ^ d. 

To evaluate the hadronic matrix elements of these operators we employ the heavy- 
quark limit mf, ^ Aqcd to get the factorization formula |6|7j 

{V^{t)mB)= \FvTI + fdidvTl\i,v)^Bmv{v)\-t (15) 

where e is the photon polarization 4-vector. Here Fy is a, B transition form factor, 
and $B, $y are leading- twist light-cone distribution amplitudes of the B meson and the 
vector meson V , respectively. These quantities are universal, nonperturbative objects. 
They describe the long-distance dynamics of the matrix elements, which is factorized 
from the perturbative, short- distance interactions expressed in the hard-scattering ker- 
nels T/ and T/^. The QCD factorization formula (fT3j) holds up to corrections of relative 
order Aqcd/"^6- 

To leading order in QCD and leading power in the heavy-quark limit, gives the 
only contribution to the B — >■ V'^ amplitude. Its matrix element is simply expressed in 
terms of the standard form factor, Tj is a purely kinematical function, and the spectator 
term T/-^ is absent. At 0{as) the operators Qi,,^ and Qs start contributing and the 
factorization formula becomes nontrivial. 

The relevant diagrams for the NLO hard-vertex corrections T/ have been computed 
in jl5ll6j to get the virtual corrections to the matrix elements for the inclusive b —>■ sj 
mode at next-to-leading order. For the exclusive modes the same corrections enter the 
perturbative type I hard-scattering kernels. The non-vanishing contributions to T/^, 
where the spectator participates in the hard scattering, are shown in Fig. ^ We can 
express both the type I and type II contributions to the matrix elements (Qi) in terms 
of the matrix element {Q7), an explicit factor of a^, and hard-scattering functions 
and Hi, which are given in |6|14j . 
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Figure 1: 0{as) contribution at leading power to the hard-scattering kernels Tf^ from 
four-quark operators Qi (left) and from Q^. The crosses indicate the places where the 
emitted photon can he attached. 

Weak annihilation contributions are suppressed by one power of AQCD/''^fe but nev- 
ertheless calculable in QCD factorization, because in the heavy-quark limit the colour- 
transparency argument applies to the emitted, highly energetic vector meson. Despite 
their suppression in Aqcd/''^65 they can be enhanced by large Wilson coefficients Ci^2 
and thus still give important corrections. This situation is relevant for B p'j. Weak 
annihilation is sensitive to the charge of the decaying B meson and thus leads to isospin- 
breaking differences between B^ p'^'y and B^ — >■ p°7. The corresponding mechanism 
is CKM suppressed in the case of i? where penguin operators give the dominant 

effect for isospin breaking. 

The total B — > V'y amplitude then can be written as 

A{B ^Vj) = ^ [A,a? + Ka^] {V^mB) (16) 

where the factorization coefficients 07(1^7) consist of the Wilson coefficient C7, the con- 
tributions from the type-I and type-II hard-scattering, and annihilation corrections. One 
finds a sizeable enhancement of the leading order value, dominated by the T^-type cor- 
rection. The net enhancement of 07 at NLO leads to a corresponding enhancement of 
the branching ratios, for fixed value of the form factor. This is illustrated in Fig. El 
where we show the residual scale dependence for B{B K*^^) and B{B~ —>■ p~7) at 
leading and next-to- leading order. As shown already in |6pi4j . our central values for 
the B —>■ K*'~f branching ratios are higher than the experimental measurements (0), 
The dominant theoretical uncertainty comes from the B —>■ Vj form factors. We used 
the light-cone sum rule (LCSR) results Fk* = 0.38 ± 0.06 and = 0.29 ± 0.04 from 
[T7] . A recent preliminary lattice QCD determination, F^* = 0.25 ± 0.05 ± 0.02 [18], 
would give a better agreement with the experimental central values. Further studies of 
heavy-to-light form factors will also benefit from developments based on factorization 
and soft-collinear effective theory (for recent discussions see |19|20|2H22p . 

Using the experimental results for the exclusive B^ K*^'y and inclusive B — > Xsj 
branching ratios together with their theory predictions, we could extract a value of the 
B K* form factor that is essentially independent of CKM factors and potential new 
physics effects. The most recent measurement of the inclusive branching ratio comes 
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Figure 2: Dependence of the branching fractions B{B'^ K*^'-)) and B{B~ — > p~7) 
on the renormahzation scale /i. The dotted hne shows the LO, the dash-dotted hne the 
NLO result including type-I corrections only and the solid line shows the complete NLO 
result. 



from the Belle collaboration, which reports 

B{B X.^r^l^^^^^ = (3.59 ± 0.47) ■ 10"^ (17) 

in the photon energy range 1.8 GeV < < 2.8 GeV. The theory predicition from a 
complete NLO QCD calculation is [16] 

BiB ^ X,7)*^',,,,«,, = (3.57 ± 0.30) ■ 10"^ (18) 

for a photon energy cutoff at Eq = 1.6 GeV. Using the approximate expression for the 
integrated branching ratio as a function of Eq in |24J, we find that 98.7% of the events 
with E^>1.Q GeV have E^ > 1.8 GeV, i.e. 

BiB ^ X,7)*^'^,,3^,, = (3.55 ± 0.30) ■ 10^^ (19) 

Since prediction and measurement are in excellent agreement for the inclusive branching 
fractions, we may consider this as a confirmation of SM short-distance physics in 6 — 57 
transitions. We can then proceed to directly extract Fk* from the measured B{B^ — >■ 
i^*07). 

With our theory prediction for the CP averaged B{B^ —>■ K*^'j) we get to very good 
approximation 

Fk* = -0.025l[!:;]?^ + 0. IbOt^ofoW 10^ B{BO ^ ir*o^)exp (20) 

Here the errors are due to the variation of the renormahzation scale mb/2 < fi < 2mb, 
which is the largest source of theoretical uncertainty [Bj. Using and adding errors in 
quadrature we get 

F^F = 0.28 ± 0.02 (21) 
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The input parameters used throughout this paper are collected in Tabled 



CKM parameters and coupling constants 


Vus 






Vub/V,k\ 




a 




0.22 


0.041 


0.09 ±0.02 


(225 ± 25) MeV 


1/137 


1.166 X lO-^GeV-2 


Parameters related to the B mesons 


rriB 


/b [25] 




rB+ 


TbO 


5.28 GeV 


(200±30) MeV 


(350 ± 150) MeV 


1.67 ps 


1.54 ps 


Parameters related to the K* meson il7j 


Fk* 






af 


af 




0.38 ±0.06 


185 MeV 


894 MeV 


0.2 ±0.2 


0.04 


218 MeV 


Parameters related to the p meson \ 


17 


Fp 


J p 


nip 


a1 




fp m 


0.29 ±0.04 


160 MeV 


770 MeV 





0.2 ±0.2 


209 MeV 


Parameters related to the uj meson 










a- 


U [26J 


0.29 ±0.04 


160 MeV 


782 MeV 





0.2 ±0.2 


187 MeV 


Quark and W-boson masses 








Mw 


(4.2 ± 0.2) GeV 


(1.3 ±0.2) GeV 


174 GeV 


80.4 GeV 



Table 1: Summary of input parameters. 



3 CKM Parameters from B{B pj)/B{B K*j) 

Ratios of different B Vj decay modes can give information on parameters in the (p, fj) 
unitarity-triangle plane with reduced hadronic uncertainties. The most natural choice is 
the ratio of the neutral B^ and B^ K*^'y branching ratios since annihilation 

effects in 5° — > p°7 are much reduced in comparison with B^ P^l- On the other 
hand, these effects can be estimated and the charged mode can also be used for a similar 
analysis. 
We define 

where the B mesons have the quark content B = {bq) and B = (hq). We will also 
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consider the CP-averaged ratios 



R 



B{B p-i) + B{B -> p7) 
B{B K*-i) + B{B K*-i) 



(23) 



Omitting the neghgible effect of direct CP violation in B ^ that is assuming 

B{B K*-f) = B{B i?*7), we may write for R = Rq, R± 



Rn 



+ R{B' 



The ratio R{B) can be expressed as 



R{B) 



td 



Vt 



ts 



t-2 



and 



«7(P7) 



R, 



R{B+)+R{B-) 



(24) 



1 -5a- 



p + irj 



p — 17] 



(25) 



Here Cp = 1/ -\/2 for p = p^ and Cp = 1 for p = p"^ 



m 



B 



m 



K* 



1.023 



and 



5a 



0-7 (pl) - o-iiPl) 



(26) 



(27) 



a^(p7) 

The coefficients in (jTTj) are understood to include the annihilation contributions. If 
annihilation effects are neglected 5a = 0{as)- The annihilation terms, on the other hand, 
contribute to 5a only at order KqcD/^b- To first approximation weak annihilation is 
induced by the leading four-quark operators and It enters the coefficients aji^pj) 
as an additive term given by |6|14j 



bu ai for B^ 
hd a2 for i?° 



p^7 
p% 



(28) 



Here ai^2 = C'1,2 + C2,i/S and 



47r^ fsfpmp 
3 FpniBmbXB 



(29) 



In the derivation of ()25|1 we have used the identity 

\n Qjij Cl>j 

Xt aS 



(30) 



and neglected the second term in the brackets in the case of -B ^ K*'~f where it amounts 
to a correction of less than 0.2% for the neutral and less than 1% for the charged mode. 
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CP averaging ()25|) and expanding in 5a we get 



R 



V, 



td 



Vt 



ts 



f-2 



«7(P7) 



1 + 2Re5a 



rf - p{l - p) 



pY + rf 



(31) 



For the case of the neutral modes (-Ro), the term proportional to Re 5a is a small cor- 
rection. The numerical value and the errors from various sources, indicated in brackets, 
are found to be 



Re 5a, = 0.002 {^^) i^s) T^^^ (fs) t^oll (F,) T^}^ («^) (32) 



The scale p has been varied between mh/2 and 2mfe and the remaining input according 
to Tabled The central value is very small because of a somewhat accidental cancellation 
between the 0{as) effects and the annihilation corrections in 5a. Adding in quadrature 
the positive and negative deviations in ()H2|1 we find 



Re5ao = 0.0 ±0.1 



(33) 



One may note further that the CKM factor multiplying Re 5a in (13111 is small for the 
region in the (p, f/) plane allowed by the standard fit of the unitarity triangle. In terms 
of the CKM angle 7 and Rh = \/ p"^ + f]"^ this factor can be written as 



f{p,v) 



V - p(i - P) 

{1-py + v 



F,2 



Rl — Rb cos 7 
— 2Rb cos 7 + -R, 



(34) 



The standard fit region, which is the most interesting for precision tests of the CKM 
framework, is roughly characterized by 



71 71 
- < 7 < - 

4 - ' - 2 



0.3 <Rb< 0.5 

This implies —0.2 < f{p,f]) < 0.2. Together with ()33p we then have 

\Re5aof{p,f])\ < 0.02 



(35) 



(36) 



This means that, under the conditions mentioned above, the correction proportional to 
Re (5a in fl31|) can be safely neglected and the relation between Rq and CKM quantities 
greatly simplifies. 

Taking into account the uncertainties from scale dependence, A^, fs, Fx*, Fp, af* 
and we get 

a^(ir*7) 

We recall that a^ is essentially free of annihilation contributions, which mainly affect a^. 
Defining 



1.01 ±0.02 



(37) 



4(P7) 
a^(fr*7) 



K 



0.98 ±0.02 



(38) 
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and recalling Q, we finally have 



Using K = 0.98, which leads to the second equality in (jH^ . this formula holds to within 
±3%. In this approximation Rt, which is the radius of a circle around the point (1, 0) in 
the {p, ff) plane, is directly given in terms of the CP-averaged ratio of branching fractions 
Rq = B{B^ p^'-j) / B{B^ — i> K*^'-)). The theoretical uncertainty is essentially reduced 
to the SU(3) breaking parameter ^ = Fk* jFp- We use the LCSR estimate ^ = 1.31 ±0.13 
[TTj . A preliminary lattice value is^ = l.l±0.1 |18j . 

In the case of the charged modes, with a decaying B^, weak annihilation dominates 
6a and we typically have 

Re5a± = -0.4±0.4 (40) 

The uncertainty is largely due to determining the strength of weak annihilation. 
This parameter is still not well known at present, but the situation can in principle be 
systematically improved |27l28j . 

The constraint in the (p, f]) plane implied by a measurement of Rq is shown in Fig. 
El For the purpose of illustration we shall assume that the results in and |1| can be 
interpreted to give 

5(5° p°7) = (0.30 ±0.12) • 10-^ (41) 

Here we have combined the average p^ / p^ /uj branching ratio from [3] and obtaining 
B{B ^ {p/uj)-i) = (0.64 ± 0.27) ■ 10"^. Dividing by 2tb+/tbo then gives (gH) as an 
estimate for B{B^ P^l)- We use the central value in (P) to compute the experimental 
ratio 

B(B^ ^ rtS) 

^o = ^(),o^^J4 = 0-007 ±0.003 (42) 

Adopting an error of ±10% for the SU(3)-breaking form factor ratio C, = 1.31 ± 0.13 
and the central value in ()42|) . we obtain the dark shaded band in Fig. El Here the 
full expression ()25|) is used, without expanding in Sa, and all theoretical parameters 
besides ^ are kept at their central values. This is justified as the theoretical uncertainty 
is entirely dominated by ^. For the same constraint, the dash-dotted lines indicate the 
la experimental uncertainty from (j42j) with fixed ^ = 1.3. 

As can be seen, the intersection of the constraints from Rq and sin 2/5 determines 
the apex (p, f/) of the unitarity triangle. For comparison, the standard fit region for 
the unitarity triangle in the (p, fj) plane [221 ^^"^ constraint from the experimental 
measurement of sin 2/5 = 0.734 ± 0.054 [201 are also shown in Fig. |21 

We finally note that the information from Rq is already becoming comparable with 
the constraint from the ratio of Bd and Bg meson mixing frequencies AM^^ and AMb^ 
[3T]. It is possible that a useful experimental measurement of Rq might actually be 
achieved before the measurement of AMb^ ■ Very interesting in this respect is the recent 
upper bound for B{B^ — >• p°7) from Babar (jH). As we have discussed above, the neutral 
mode is favoured theoretically because of the small impact of annihilation effects. In 
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0.5 



1 



P 



Figure 3: Constraints implied by Rq in the (p, r]) plane. The experimental value used 
is Rq = 0.007 ± 0.003 (see text for further explanation). The width of the dark band 
reflects a ±10% variation of ^ for central Ro- The dash-dotted lines display the error 
from Ro while ^ is kept at its central value. The region obtained from a standard fit 
of the unitarity triangle (irregularly shaped area) and the constraint from sin 2/? (light 
shaded band) are overlaid. 

addition, it turns out that the upper bound for B{B^ —>■ p^'j) is particularly strong in 
comparison with the bound for the charged mode (0), even after correcting for an isospin 
factor of 2 and the B^/B^ lifetime difference. We thus prefer to use the neutral mode 
directly for placing an upper bound on Rt, rather than the combined result of the three 
modes (p^/p°/cc;)7, which was the choice made in [Hj. Using ()39j) . the recent Babar limit 
dH) together with (jT} implies 



Rt < 0.82 — 
1.3 



(43) 



This is equivalent to 



Vtd 




Vtd\ < 7.3 ■ 10"^ — 
' 1.3 



(44) 



Vts 
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0.5 1 

P 



Figure 4: Upper bound on Rt (the distance from point (1,0)) implied by the Babar 
hmit B{B^ p^^) < 0.4 ■ 10^® in the {p,fi) plane. The curves correspond (from left to 
right) to ^ = Fx*/Fp = 1.4, 1.3 and 1.2. The region obtained from a standard fit of the 
unitarity triangle (irregularly shaped area) and the constraint from sin 2(3 (light shaded 
band) are overlaid. 

The bound may be compared with the 2a range 

6.5-10-^ <\Vtd\< 9.5 -10-^ (45) 

obtained from a standard fit of the unitarity triangle jSHj- For 30% SU(3) breaking 
in the ratio of form factors, ^ = 1.3, more than half of the range ()45|) is excluded by 
(jUI). Should the amount of SU(3) breaking be less than 30%, the bound would be even 
stronger. An illustration of the Babar bound in the {p, ff) plane is given in Fig. |3] 

4 B ^ and B plu 

As we have seen, the rare decay i? — > p7 is a clean probe of flavour physics, except 
for the sizable uncertainty in the form factor Fp. Other uncertainties are quite well 
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under control within a treatment of the decay at next-to-leading order in QCD and a 
leading-order evaluation of power-suppressed annihilation effects. This is the case in 
particular for the neutral channel 5° where weak-annihilation effects are small. 

The sensitivity to Fp can be reduced by taking the ratio of 5 — > p7 and B K*^ 
branching fractions, as we have discussed in the previous section. Then the impact of 
long-distance hadronic physics is limited to SU(3) breaking in the ratio C, = Fx*/Fp. 
While this is certainly an advantage, the exact deviation from the SU(3) limit = 1 
remains at present a significant source of uncertainty. 

In this section we discuss a possibility to reduce hadronic uncertainties in a different 
way, using the ratio of -B — plv and B ^ p'y decay rates. The simplification occurs 
because relations exist between the corresponding form factors in the large energy limit. 
Since only B ^ p transitions are involved, the problems with SU(3) breaking are avoided 
and only isospin symmetry needs to be assumed, which should be valid to within a few 
percent. The existence of relations between the form factors in the large energy limit and 
their potential usefulness for phenomenology were first pointed out in The results 
of [32j were put on a field theoretical basis within the soft-coUinear effective theory 
(SCET) and extended to higher order in QCD |33l22l2lj . These relations were applied 
to extract information on the form factor in i? ^ K*'~f for use in other channels such 
as B ^ K*l^l~ [34J. Previously, the authors of have investigated the possibility to 
relate B plv in a certain region of phase space with B — > K*'^. This suggestion is 
similar in spirit to our proposal, but the analysis of |^ was based only on the heavy 
quark limit, instead of the full large energy relations from and was still affected by 
SU(3) breaking. In addition, our discussion also includes short-distance QCD corrections 
at next-to-leading order. 

The differential decay rate for B ^ plu is given by 



d^T{B plv) 
ds dz 



B 



P 2567r3 



'l-z) 



-Hi + 



+ 



(46) 



Here 



and 



w = w{s, r) 



1 + + - 2s - 2r - 2sr 



(47) 



z = cose s = ^ r = —f^ 0.021 (48) 

where is the dilepton invariant mass and 9 is the angle between the momenta of the 
neutrino and the B meson in the dilepton centre-of-mass frame. Equivalently, 9 is the 
angle between the charged-lepton momentum and the direction anti-parallel to the B 
momentum in this frame. The same definition of 9 is valid for B — »• plv with either a 
positive or a negative charged lepton. The kinematical range for s and z is 



< s < (1 - v/r)^ 



- 1< z < 1 



(49) 



The Hi = Hi{s) are helicity form factors. They can be expressed in terms of the vector 
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and axial vector form factors V{s), Ai{s) and A2{s) as 



Ha 



1 



V 



2v^ 



1 + v/r 
:i-s-r)(l + v^)Ai 



w 



1 + 



(50) 
(51) 



where we use the conventions of |321 for V, Ai, A2, which, in particular, are positive real 
quantities. 

The CP-averaged decay rate for B p'-f can be written as 



T{B ^ p7) 



G^am^ml 



327r^ 



(52) 



where we have used the approximation, explained in the previous section, that corre- 
sponds to neglecting the term ~ Re (5a in pip . As we have seen, this is a very good 
approximation for 5° — p°7. If a more accurate treatment is desired, or the analy- 
sis should be applied to B"^ p'^'j, the following discussion can be generalized in a 
straightforward way using the complete expression based on (fTB|). Combining and 
(|52|) we find for the case of neutral B mesons 



k{s, z) 



VudYub 



VtdVtb 



Aa ml (1 -r)3|KdP|a7(p7)P dB{B^ ^ plu) 



TT ni l 



B{BO ^ PO7) 



ds dz 



(53) 



Here B^ plv can be either one of the two channels B^ —>■ p^l v 01 B^ ^ p l^v and / 
may be an electron or a muon. The hadronic quantity k is defined as 



k{s, z) = w{s, r)^/^ 



(l-^r Hj (1 + HI Hi 
2 Py 2 F2 ^ ^ ' Fl 



(54) 



The differential branching ratio and the function k in ()53|) may be replaced by their 
integrated versions 



AB{s,e) 
where 



ds 



^ , dB(B^ plv) 
dz- 

l-e 



ds dz 



< s < 



K{s, 



ds 



dzk{s,z) (55) 



= (1 - Vrf ^ 0.73 < e < 2 (56) 

such that the fully integrated branching fraction for 5*^ plv is given by AB{smax, 2). 
The relation ()53j) allows us to determine the CKM parameter 



yudVub 



VuVth 



;i - pf + r] 



F,2 



(57) 



in terms of observable B p'j and B plv branching fractions, known quantities, 
and the hadronic function k. The main virtue of this expression is that in the large 
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energy limit hadronic form factors cancel in the ratios H±/Fp. This is not the case 
for Ho/Fp, but its contribution can be suppressed by selecting events in the vicinity of 
z = 1. As a consequence, ()5Hj] can be turned into a theoretically clean expression for the 
determination of the CKM ratio in (jHTj) . 

In the large energy limit the form factors H±, Hq and Fp can be written in terms of 
just two independent form factors C±{s) and ClK-^) using ^32j 



A,{s) 



1 — s + r 
1 + 



C±{s) V{s) = {l + V^)a{s) 



A2{s) = (1 + v^) 



C±is) 



2v^ 



Cii(^) 



1 — s + r 

Together with (|50|l . (j5T| these relations imply 

H±{s) = s + vtV^] C±{s) 



F,^ri(o) = a(o) 



Ho{s) = V^{l + s-r)C±{s) + 



w 



s + r 



Cllis) 



(58) 

(59) 

(60) 
(61) 



These results are valid in the heavy-quark limit and the limit of large energy of the 
recoiling p-meson 



me 
2 



s + r) 



(62) 



In this approximation the ratio / is independent of hadronic form factors. For not 
too large values of s this ratio is strongly suppressed, Hj^/ = 0{r). More importantly, 
also H_{s) and Fp depend on the same form factor C±('5); which is to be evaluated at 
s = in the latter case. As a consequence, we may write 



4s(l + 6is + ...) 



(63) 



expanding the form factor ratio in a Taylor series. The leading term in this ratio for 
small s is largely free of hadronic uncertainties in the large energy limit. The higher- 
order corrections only depend on the shape of C±("5)5 not on its absolute normalization, 
and can in principle be determined from a fit to the shape of the observed spectrum in 
s. The coefficient bi is related to the slope of (± and can be written as 



Ci(0) 

a(o) 



1 



1 — r 



(64) 



When fitting the ratio in ()6Hj] to the experimental spectrum, other parametrizations for 
the shape may, of course, be chosen. The Taylor series 1 + bis + . . . could be replaced 
for instance by the pole form 1/(1 — /9is)^, or a combination of the two. 

In [33j the corrections of order have been computed to the relations between form 
factors in the large energy limit. There is no relative correction between Ai and V to all 
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orders in Therefore the correction of the ratio V/Fp given in fSBj also apphes to 

H_/Fp. Taking these effects into account, the leading term (4s) in is modified to 



/ 2a,(/ii) 


"l + 21ni^" 


V 37r 


m6_ 



37r 1/(0) 



:i + v^)^ (65) 



where the first term with Osifii) ~ 0.22 refers to the vertex correction and the second 
with as(/i2) ~ 0.34 to the hard spectator interaction. The usual renormalization scheme 
of the form factor Fp, adopted in this paper and used in (jHSI), corresponds to the MS 
scheme with anticommuting 75 (NDR). Numerically, the QCD correction factor amounts 
to (1 — (0.15 ± 0.10)) using the estimates in |^. The dominant uncertainty comes 
from AF_i_, which depends on properties of the 5-meson light-cone wave function. This 
quantity is poorly known at present, but improvements should be possible in the future 
and would lead to a reduction in the uncertainty. 

It is interesting to compare the above analysis with the results for the form factors 
obtained using the method of light-cone QCD sum rules ^Zj. With the form factors 
computed in J7j one finds 



H 



2 I 



F2 

p 



4.25s (1 + 0. 59s + 0.65s^ + ...) (66) 



The leading term agrees very well with the prediction at leading order in the large energy 
limit (j63p . Taking the QCD corrections into account according to (j65|) . the prediction 
for this term is typically about 15% lower. On the other hand, the result in (jUBj) has an 
uncertainty of about 30% |T7]. Nevertheless, the general level of agreement of the sum 
rule calculations, which include subleading corrections in l/rrif,, with the large energy 
limit, is consistent with the assumption that power corrections are of moderate size. 

In contrast to H±, the longitudinal form factor Hq is dominated by Cy, which is not 
cancelled in the ratio Ho/Fp. The third term in (jS^ can still be estimated theoretically, 
but will be affected by larger uncertainites. As mentioned above, in order to reduce 
its importance a cut on the angular variable z may be imposed, restricting z to be in 
the vicinity of +1 or —1. The latter case is not interesting, since it would strongly 
suppress the if_ contribution, leaving only the contribution from which is very 
small. Parametrizing the cut below z = 1 by e as defined in ()55|) and performing the 
angular integration, we find for K{s, e) 



K{s, e) = / ds \fw 
Jo 



6 F2 V Q J Fp \ 3 7 F2 



(67) 



The full angular range is obtained for e = 2 and in this case all three e-dependent 
coefficients become equal to 4/3. For small e, on the other hand, a strong hierarchy 
exists, which is clearly visible in (jU7j) . The contribution from Hq is suppressed with 
respect to the if_-term by a factor of e/2, that is by one order of magnitude for e = 0.2. 
The corresponding suppression of the if_|_-term is even by a factor of 6^/12, in addition 
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to the fact that H^/H_ is aheady small for moderate values of s. The contribution from 
iJ+ is therefore entirely negligible in the following discussion. 
Neglecting all terms of C(e^), (jHTj) simplifies to 



K{s, e) = / ds \fw 
Jo 



ZJ2 tt2 -\ 

n r, rir, — n 



p2 

P P 



(68) 



The validity of the large energy limit, with the model- independent normalization of 
H^/Fp in (jUHj) . ()65|) . requires moderate values of s. Enhancing this term in ()67p requires 
small e. As a typical example one may concentrate on the part of phase space defined by 
s = 0.4 and e = 0.2. The relative number of — > plu events in this region (0 < s < 0.4, 
0.8 < -2 < 1) is given by 

ir(0.4,0.2) , , 

For this estimate we have evaluated K{s,e) employing the form factors from 17j. A 
measurement of B ^ plu has been reported by CLEO ^37) . 



5(5° ^ pZz/) = (2.17 ±0.73) ■ 10"^ (70) 

and BaBar ^ : 

B{B° plu) = (2.57 ± 0.79) ■ 10"^ (71) 

The effective branching ratio of i? — >^ pip events in the above region of phase space would 
then be about 10~^. 

The first term in (|68|) is determined by the measured shape of the s-distribution and 
the mo del- independent normalization in (j63p . (j65|) . The small correction from the second 
term in ()68p could either be estimated theoretically, or be isolated in the data by varying 
e. With the form factors from J7] we have for instance 

ii'(0.4,e) = 0.57 e + 0.25 (72) 

Once K{s, e) is known, the measured values of Ai?(s, e) and B{B^ p°7) determine 
the CKM quantity in ()57|) using ()53|) . This CKM ratio provides us with an interesting 
constraint in the (p, f/) plane, which is illustrated in Fig. for a hypothetical measure- 
ment of iVudVub/Vidl"^ = 0.16 ± 0.04. We observe that the constraint is quite stringent, 
in particular in the important region corresponding to the standard fit results, and even 
for the rather moderate precision of ±25%. 



5 Isospin Breaking in B 

The CP averaged isospin breaking ratio can be defined as 



A(r,) = nB^^v^i)-vnB-^v-,) 



16 



0.6 



0.4 
0.2 
isr 
-0.2 
-0.4 



-0.6 

-0.8 -0.6 -0.4 -0.2 0.2 0.4 

P 

Figure 5: Constraints in the (p, r/) plane implied by iKdKtfe/VjdP = 0.16 ± 0.04. For 
comparison, the standard fit region is indicated by the shaded area. 

with w = 1 for = K* and v = 1/2 for V = p. This ratio has a reduced sensitivity to 
the nonperturbative form factors. As already discussed, in our approximations, isospin 
breaking is generated by weak annihilation contributions. Kagan and Neubert found a 
large effect from the penguin operator Qq on the isospin asymmetry A{K*'^) jSH]. Our 
prediction /S.[K*'^) = (3.91*^^9)% (see flj) is in agreement with the experimental results 
(Belle in 0, Babar 40 ) 

A{K*-f) = +0.034 ± 0.044 ± 0.026 ± 0.025 (Belle) (74) 
A{K*-f) = +0.051 ± 0.044 ± 0.023 ± 0.024 (Babar) (75) 

Here the errors are statistical, systematic and from the B^/B^ production ratio. 

For S ^ p7 we find a strong dependence of the isospin asymmetry on the angle 7 of 
the unitarity triangle. As seen in Fig. IHl the 7 dependence is in particular pronounced 
for the zero crossing of A(p7) around 7 = 60°, the value favoured by the standard UT 
fits. 

Once a measurement of both the charged and neutral B —>■ pj modes is available, 
the isospin-asymmetry A(p7) can be used to constrain the unitarity triangle. For the 
purpose of illustration we plot in Fig. [7[ in addition to the Rq and sin 2/3 bands shown 
already in Fig. |21 the implication of an assumed measurement of A(p7)e^p = 0, which 
would correspond to the Standard Model prediction for a CKM angle 7 = 60°. The 
dominant theoretical uncertainty comes from the hadronic parameter Xb and from the 
variation of the renormalization scale. 
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Figure 6: The isospin-breaking asymmetries A(p7) and A{K*'~f) as a function of the 
CKM angle 7 with (sohd) and without (dashed) the inclusion of QCD penguin operator 
effects. 

6 B ^ (J J 

In this section we briefly consider the decay u^'y and discuss differences to the 

related mode 5° p°7. We consider p° and uj^ as pure isospin-1 and isqspin-0 modes, 
respectively, and neglect p — u mixing. We use the convention oj^ = = ^^^J^ ■ 

To leading order in the heavy-quark limit and next-to-leading order in «<, both the 
p° and uj^ meson in i? ^ 1^7 are produced from a dd pair. Therefore, to get the 
B —y uj^'j decay amplitude, we can use the one for B — >■ p^j with obvious replacements 
for the vector meson decay constant, mass, LCDA and form factor in the factorization 
coefficients a^^^{p^'^) |6ll2ll4j . The relevant input parameters for all the decay modes 
are compiled in Tabled 

A few comments are in order. The best known input parameter for the vector mesons 
is the mass, which can be found in the Review of Particle Physics [41 . Using r-decay data 
and the purely leptonic decay modes of p° and uj^ one can extract the respective decay 
constants /p with negligible uncertainty [21]. The other vector meson parameters, such 
as the B ^ V form factors were taken from QCD sum rule estimates ^Tj. We take the 
same values for the u and p mesons, which should be a reasonable assumption, even 
though this equality could be broken by Zweig-rule violating effects. The latter are, 
however, suppressed by l/Nc. For instance, the decay constants fp and /^^ differ by 10%. 

In |6ll4j we included weak annihilation contributions to i? ^ Vj although they are 
suppressed by one power of AQcn/^b- The reason for including these power-suppressed 
contributions was that they are in part enhanced by large Wilson coefficients, they are 
calculable in QCD factorization and they can be used to estimate isospin-breaking effects. 
For B — >■ u'j annihilation contributions are also calculable and they are the source of 
specific differences (apart from form factors) between B^ p°7 and B^ ui^'-f. Those 
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Figure 7: Same as Fig. 01 including the implication of a measurement of A{frf)exp = 
(curved band on the right). The width of the band reflects the theoretical uncertainties 
from varying the hadronic parameter and the renormalization scale /i. (The effect of 
isospin breaking in the form factors is neglected here.) 



are due to the fact that p° and are isospin-0 and isospin- 1 states, respectively. In 
the following we will use the notation of section 4.5 in If, in figure |H1 the photon 
emission is from the light quark in the B meson, the annihilation amplitude contains 

^ 27r^ fBmyfv 
Fy niBmbXB 

whereas the (55,6 insertion with the photon emitted from one of the vector meson con- 
stituent quarks leads to 



V ty mBTnf, Jq i.^) 



The annihilation coefficients for uj^'j then are 

Cn(^°7) = Qd [+a2&" - 2b-{a, + a,) + a,b^ + + rf^)] (78) 

<nn(^°7) = Qd [-26"(a4 + ae) + a,b^ + a^id^ + f/^)] (79) 
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Figure 8: Annihilation contribution to the B — > V'-f decay. The dominant mechanism 
for operators Qi,,a is the radiation of the photon from the hght quark in the B meson, as 
shown. This amphtude is suppressed by one power of lyqcD/^b-i but it is still calculable 
in QCD factorization. Radiation of the photon from the remaining three quark lines is 
suppressed by {^QCD/''^hf' for operators Qi.,a- For operators Qs^e, however, radiation 
from the final state quarks is again of order h.qcol'^h- 

The difference compared to aiJnni^p^^^ is the sign change of the 02 contribution and the 
additional isospin-0 contribution —26'^ (04 + oe). 

Numerically the lo and (P annihilation coefficients are given by 



1 



+0.0268 +0.0281 -0.0060 +0.0446 
+026"^ -26'^ (04 + 06) +046^ +a6« + ci^) 

= -0.0028 



= -0.0312 (80) 



-0.0296 -0.0066 +0.0446 
-026'' ^a^V +a6« + (i^) 

For comparison we quote the corresponding numbers for the p~7 channel 

2 



(81) 



^{P 7) = 



+0.2902 -0.0066 -0.0112 +0.0223] = +0.1965 

+ai6^ +046'' ^QslQua%d% ^a^d% 



(82) 



where the annihilation component is considerably larger. This has to be compared with 
a^{oj^^ — ai^{p^^) — —0.4154 — 0.0685i, to which the annihilation coefficients are added. 
For central values of all input parameters, fi = mi,, and our default choice for the CKM 
angle 7 = 58°, we get the following CP-averaged branching ratios: 

(83) 
(84) 

(85) 

Within the parametric and theoretical uncertainties the B^ wf and B^ p7 branch- 
ing ratios can be considered equal, neglecting any possile difference in the respective 
form factors. 



S(S° - 




= 0.84 


10 


B{B°- 


-A) 


= 0.81 


10 


B{B^ - 


.pS) 


= 1.81 


10 
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7 Conclusions and Outlook 



We have studied constraints on the CKM unitarity triangle from observables in the 
exclusive radiative decays B K*'y, B — >■ p7, and B wy, as well as the exclusive 
semileptonic decay B plv. Within the framework of QCD factorization we have 
worked at next-to-leading order in to leading order in the heavy-quark limit. Power 
corrections from weak annihilation have also been included. Important information on 
the unitarity-triangle parameters p and f/ can be obtained from the ratio Rq of the 
neutral B^ —>■ p^j and B^ — >■ K*^^ branching ratios. This ratio measures to very good 
approximation the side Rt of the standard unitarity triangle. Annihilation effects are 
negligible in this case. The theoretical uncertainty in the relation to Rt comes in essence 
solely from the form-factor ratio C, = Fx*/Fp, which differs from unity only because of 
SU(3)-breaking effects. Using the latest bound on B{B^ — ^ from Babar we find 
Rt < 0.81 (C/1.3) or \Vtd\ < 7.3 ■ lO'^ (^/1.3) (see also Fig. ^. 

Similar constraints in the (p, ff) plane can be obtained from the isospin asymmetry 
A(p7) once a measurement of this quantity is available. 

We propose to gain complementary information in the (p, ff) plane through the 
B plv and B ^ p'y decay rates, which can be related to the CKM parameter 
For events where the momenta of the neutrino and the B meson 
are parallel in the dilepton centre-of-mass frame, this relation is free of hadronic form 
factors in the large energy limit. This allows a theoretically clean determination of the 
above CKM ratio. We have shown that even a moderate experimental precision can 
yield a stringent constraint in the (p, f/) plane. 

Finally, we have calculated the annihilation effects in the B wy decay amplitude 
which turn out to be very small. 

An improved determination of the B ^ V form factors and, in particular, the form- 
factor ratio ^, remains an important task for the future. More precise experimental 
measurements, specifically the individual measurements of B{B^ p^'-f) and B{B~^ 
p'^'j) are eagerly awaited. These measurements can lead to results on Rt competitive 
with those from Bg-Bg mixing. An experimental analysis of the differential B — ^ plu to 
B ^ pj decay rate ratio can circumvent the form-factor related uncertainties to a large 
extent and will thus be of particular interest. 
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of the B —>■ uj form factor. This research was supported in part by the National Science 
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